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Abstract 

In  this  paper  we  discuss  controllability  of  Lie-Poisson  reduced  dynamics  of  a 
class  of  mechanical  systems.  We  prove  conditions  (boundedness  of  coadjoint  orbits 
and  existence  of  a  radially  unbounded  Lyapunov  function)  under  which  the  drift 
vector  field  (of  the  reduced  system)  is  weakly  positively  Poisson  stable  (WPPS). 
The  WPPS  nature  of  the  drift  vector  field  along  with  the  Lie  algebra  rank  condition 
is  used  to  show  controllability  of  the  reduced  system.  We  discuss  the  dynamics, 
Lie-Poisson  reduction,  and  controllability  of  hovercraft,  spacecraft  and  underwater 
vehicles,  all  treated  as  rigid  bodies. 


1  Introduction 


A  geometric  approach  to  the  study  of  mechanical  systems  has  had  a  profound  influ¬ 
ence  in  recent  years  on  our  understanding  of  dynamics  and  control  aspects.  Playing 
an  essential  role  in  this  are  recent  developments  in  reduction  theory  (cf.  [1]),  i.e.  the 
exploitation  of  invariance  of  the  controlled  dynamics  to  a  group  of  transformations  (the 
symmetry  group).  The  existence  of  a  symmetry  group  permits  dropping  of  the  dynamics 
to  a  lower  dimensional  (reduced)  space.  Lagrangian  reduction  [2,  3]  involves  dropping 
the  Euler-Lagrange  equations  to  the  quotient  of  the  velocity  phase  space  given  by  the 
symmetry  group  while  Hamiltonian  reduction  involves  projecting  the  Poisson  bracket  to 
the  reduced  (quotient)  space  which  also  inherits  a  Poisson  structure  (attributed  to  Lie 
and  Berezin-Ivirillov-Kostant-Souriau,  see  the  work  of  Weinstein  for  historical  remarks 
[4]).  In  particular,  if  the  configuration  space  of  the  system  can  be  identified  with  a  Lie 
group  G ,  a  left  invariant  Hamiltonian  on  T*G  gives  rise  to  reduced  dynamics  on  T*G/G 

*This  research  was  supported  in  parts  by  grants  from  the  National  Science  Foundation’s  Engineer¬ 
ing  Research  Centers  Program:  NSFD  CDR  8803012  and  by  the  Army  Research  Office  under  Smart 
Structures  URI  Contract  No.  DAAL03-92-G0121. 
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which  is  isomorphic  to  Q*  the  dual  of  the  Lie  algebra  of  G.  The  reduced  bracket  is  now 
the  Lie-Poisson  bracket.  The  complete  dynamics  can  then  be  reconstructed  from  the 
reduced  system.  The  geometric  phase  [5]  associated  with  a  trajectory  of  the  reduced 
system  describes  the  motion  of  the  complete  (lifted)  system.  There  has  also  been  recent 
progress  in  the  area  of  control  in  the  presence  of  nonholonomic  constraints  [6,  7,  3,  8]. 

For  a  large  class  of  mechanical  systems,  the  configuration  space  can  be  identified  with 
a  Lie  group  G.  Often  the  dynamics  of  such  systems  are  G- invariant  and  hence  they 
can  be  reduced  to  obtain  a  set  of  Lie-Poisson  reduced  dynamics  on  T*G/G.  Examples 
of  such  systems  include  hovercraft,  spacecraft  and  underwater  vehicles  modeled  as  rigid 
bodies.  The  design  and  control  of  autonomous  versions  of  these  vehicles  has  been  of 
much  recent  interest.  For  example  the  amphibious  versatility  of  hovercraft  has  given 
them  a  role  in  specialized  applications  including  search  and  rescue,  emergency  medical 
services,  ice  breaking,  Arctic  off-shore  exploration,  and  recreational  activities  [9].  Certain 
environmental  aspects  (such  as  ice-roughness,  Arctic  rubble  fields  etc.)  also  provide  a 
niche  for  operations  by  hovercraft.  Similarly  a  growing  industry  in  underwater  vehicles 
for  deep  sea  explorations  has  led  to  the  demand  for  more  versatile,  robust  and  high 
performance  autonomous  vehicles  that  can  cope  with  actuator  failures,  disturbances, 
exploit  sensor  based  local  navigation  etc. 

In  this  paper  we  discuss  the  controllability  of  the  Lie-Poisson  reduced  dynamics  of  a 
class  of  mechanical  systems  which  include  as  examples  hovercraft,  spacecraft  and  bottom 
heavy  underwater  vehicles.  In  each  case  we  identify  the  configuration  space  with  a  Lie 
group  G.  The  G-invariance  of  the  Hamiltonian  and  the  forcing  term  (control)  is  used  to 
obtain  a  set  of  Lie-Poisson  equations  on  T*G/G  which  is  isomorphic  to  Q *  the  dual  of 
the  Lie  algebra  of  G.  We  show  that  depending  on  the  existence  of  a  radially  unbounded 
Lyapunov  type  function,  the  drift  vector  field  (of  the  reduced  system)  is  weakly  positively 
Poisson  stable  (WPPS).  The  WPPS  nature  of  the  drift  vector  field  along  with  the  Lie 
algebra  rank  condition  is  used  to  show  controllability  of  the  reduced  system. 

The  paper  is  organized  as  follows.  In  section  2  we  present  a  brief  overview  of  Lie-Poisson 
reduction.  In  section  3  we  present  our  main  result  on  controllability  of  Lie-Poisson 
reduced  dynamics.  In  section  4  we  discuss  in  some  detail,  the  dynamics,  reduction  and 
reduced  space  controllability  of  the  hovercraft  ,  the  spacecraft  and  the  underwater  vehicle. 
Conclusions  and  future  work  is  discussed  in  section  5. 


2  Lie-Poisson  Reduction 


Recall  (cf.  [1,  10] )  that  if  G  is  a  symmetry  group  acting  on  a  Poisson  manifold  M,  then  the 
quotient  manifold  M/G  inherits  a  Poisson  structure  so  that  whenever  P,  Q  :  M/G  — >■  R 
correspond  to  G  invariant  functions,  P,  Q  :  M  — >■  R,  their  Poisson  bracket  {P,  Q}m/g 
corresponds  to  the  G-invariant  function  { P  Q }  a /  •  If  H  :  M  — >■  R  is  a  G-invariant 
Hamiltonian,  then  H  descends  to  H  :  M/G  — >■  R  and  determines  the  reduced  dynamics 
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on  M/G.  The  solutions  of  the  reduced  Hamiltonian  system  on  M/G  are  projections  of 

the  solutions  of  the  complete  system  defined  on  M.  In  particular  if  M  =  T*G  and 
G  acts  on  itself  by  left  translations  then  M/G  =  Q*,  the  dual  of  the  Lie  algebra  of  G, 
a  left  invariant  Hamiltonian  on  T*G  gives  rise  to  reduced  dynamics  on  Q'_  (the  space 
Q *  associated  with  the  minus  Poisson  structure).  The  reduced  bracket  is  now  the  minus 
Lie-Poisson  bracket,  defined  in  its  coordinate  free  from  by 


{F,H}_{. r)  =  -  <x,[VF(x),VH(x)]  >  . 

Let  { A" | .  •  •  • ,  A,.}  and  { A" \ ,  •  •  • ,  X/ }  be  a  basis  for  the  Lie  algebra  Q  and  the  dual  basis 

Q*  respectively,  i.e.  <  Xf,Xj  >=  5*.  Any  //,  G  Q*_  can  be  expressed  as  ^ 

and  the  Lie-Poisson  bracket  of  two  differentiable  functions  P,  Q  G  C00(^*)  is  given  by 


{F,P}_(h) 


f  ckuk^d— 
h=i 13  °ii:'  °ij] 


(i) 


where  c*.,  j,  k  =  1,  •  •  • ,  r  are  the  structure  constants  of  Q  relative  to  a  basis  {Ad,  •  •  • ,  Ad). 
Equivalently  (1)  can  be  written  as 

{F1H}_(ri  =  \7FTA(v)\7H  (2) 


where 

'.Mli/ij  E  4/ 
k= 1 


The  rank  of  Poisson  tensor  A  determines  the  nontrivial  Casimirs  of  //.  The  Lie-Poisson 
reduced  dynamics  can  now  be  expressed 

fii  =  i  =  1,  ■  ■  ■  ,r  (3) 

where  H  is  the  reduced  Hamiltonian. 

In  the  rest  of  this  paper  we  will  denote  the  Lie-Poisson  reduced  dynamics 
with  the  following  notation. 


//  =  {//.P}-  (4) 

The  /tli  component  on  the  right  hand  side  being  {//,. ,  Lf}_  ,  /  =  1,  •  •  • ,  r. 

The  induced  symplectic  foliation  by  Lie-Poisson  bracket  on  Q*  has  a  particularly  nice 
interpretation  in  terms  of  the  dual  to  the  adjoint  representation  of  the  underlying  Lie 
group  G  on  the  Lie  algebra  Q.  This  is  given  by  the  following  theorem,  which  appears 
to  be  due  to  Kirillov[ll,  12],  Arnold  [13],  Kostant  [14]  and  Souriau  [15],  though  similar 
ideas  first  appear  in  the  work  of  Lie,  Borel  and  Weil. 
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Theorem  1  Let  G  be  a  connected  Lie  group  with  coadjoint  representation  Ad*G  on  Q* . 
Then  the  orbits  of  Ad*G  are  immersed  submanifolds  of  Q *  and  are  precisely  the  leaves 
of  the  symplectic  foliation  induced  by  the  Lie-Poisson  bracket  on  Q* .  Moreover,  for  each 
g  G  G,  the  coadjoint  map  Ad*G  is  a  Poisson  mapping  on  Q *  preserving  the  leaves  of  the 
foliation. 


As  flows  of  (4)  remain  on  coadjoint  orbits  on  which  they  started,  the  geometry  of  coadjoint 
orbits  plays  an  important  role  in  understanding  the  dynamics  of  the  Lie-Poisson  reduced 
equations. 


3  Poisson  Stability  and  Controllability 


The  state  space  of  a  large  class  of  mechanical  systems  such  as  hovercraft,  spacecraft 
underwater  vehicle  etc.  can  be  identified  with  a  Lie  group  G.  The  Hamiltonian  dynamics 
(defined  on  T*G )  of  these  systems  subject  to  external  forces  can  be  written  in  the  form 
of  a  control  system  as 


m 

X  =  f{x)  +  Y^9i(x)ui  (5) 

i—  1 

where  x  G  T*G ,  f(x)  =  {x,H}  and  u  =  (rq,  •  •  -um).  ( H  is  the  Hamiltonian  defined  on 
T*G ).  Here  we  do  not  assume  that  g,  are  Hamiltonian  vector  fields.  Often  we  observe 
that  the  vector  fields  /  and  gf  s  are  G- invariant.  This  allows  us  to  drop  the  the  vector 
fields  /  and  g,  from  T*G  to  T*G/G  and  the  reduced  dynamics  take  the  form 


m 

h  =  fih)  +  J2di(h)ui  (6) 

i—  1 

where  //,  G  T*G/G ,  /  and  g,  are  the  projections  of  /  and  g  on  T*G/G.  From  the  discussion 
in  section  2  we  know  that  /  =  {fi,  iL}_  where  H  is  the  reduced  Hamiltonian.  Hence  (6) 
can  be  written  as 


m 

/'  =  UlH}-  + 

i—  1 


(7) 


Studying  controllability  of  systems  of  the  form  (7)  or  of  more  general  systems  of  the  form 


m 

x  =  f(x)  +  Hi  (-All;,  x  Gl"  u=  («1,  •  •  •  um)  G  U 

i—  1 


(8) 
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is  usually  a  hard  problem.  We  know  that  if  a  system  of  the  form  (8)  satisfies  the  Lie 
algebra  rank  condition  (LARC)  then  it  is  locally  accessible,  and  in  addition  if  /  =  0 
then  LARC  implies  that  the  system  is  controllable.  (See  appendix  for  details.)  While 
the  kinematic  equations  of  motion  can  often  be  written  as  a  drift  free  system,  once 
dynamics  are  included  LARC  does  not  imply  controllability.  Proving  controllability  is 
usually  much  harder  than  proving  accessibility.  In  [16]  sufficient  conditions  are  given,  in 
terms  of  a  “group  action”,  that  a  locally  accessible  system  is  also  locally  reachable.  In 
[17]  sufficient  conditions  for  the  controllability  of  a  conservative  dynamical  polysystem  on 
a  compact  Riemannian  manifold  are  presented.  More  recently  this  result  was  extended 
by  [18]  to  dynamical  polysystems  where  the  drift  vector  field  was  required  to  be  weakly 
positively  Poisson  stable.  We  extend  this  result  to  Lie-Poisson  reduced  dynamics.  We 
prove  conditions  under  which  the  reduced  dynamics  are  controllable.  Before  we  present 
our  results  we  introduce  some  definitions  and  related  theorems  regarding  Poisson  stable 
systems.  We  follow  the  development  in  [18,  19,  20,  21,  22,  23]. 

Let  A"  be  a  smooth  complete  vector  field  on  M  and  let  <pf  (•)  denote  its  flow. 

Definition:  A  point  p  G  M  is  called  positively  Poisson  stable  for  X  if  for  all  T  >  0  and 
any  neighborhood  Vp  of  p,  there  exists  a  time  £  >  T,  such  that  <pf{p)  G  Vp.  The  vector 
field  X  is  called  positively  Poisson  stable  if  the  set  of  Poisson  stable  points  for  X  is  dense 
in  M. 

Definition:  A  point  p  G  M  is  called  nonwandering  point  of  X  if  for  all  T  >  0,  any 
neighborhood  Vp  of  p,  there  exists  a  time  t.  >  T  such  that  4>f  (Vp)  f]Vp  ^  <p,  where 
tf  (Vv)  =  {p(q)\  qeVv}. 

One  should  observe  here  that  though  it  is  a  sufficient  condition  that  the  nonwandering 
set  of  a  positively  Poisson  stable  vector  field  is  the  entire  manifold  M,  there  could  exist 
weaker  conditions  under  which  the  nonwandering  set  is  M.  This  gives  rise  to  the  definition 
of  a  weakly  positively  Poisson  stable  (WPPS). 

Definition:  The  vector  field  A"  is  called  weakly  positively  Poisson  stable  if  its  nonwan¬ 
dering  set  is  M. 

The  following  theorem  on  controllability  is  due  to  Kuang-Yow  Lian  et.  all.  [18].  Earlier 
versions  of  this  theorem  and  the  corollary  that  follows,  where  the  hypothesis  required  / 
to  be  only  Poisson  Stable,  are  due  to  Lobrv  [17],  Bonnard  and  Crouch  [24], 

Theorem  2  If  the  system 


m 

X  =  f(x )  +  gi(x)ui ,  u  =  (uu  ■  ■  ■  Um)  G  U  C  HP" 

i—  1 

ivhere  U  contains  {u  \  \uf\  <  Mi  ^  0,  /,•••,  m]  is  such  that  f  is  a  weakly  positively 
Poisson  stable  vector  field,  then  the  system  is  controllable  if  the  accessibility  LARC  is 
satisfied. 
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Corollary  1  If  the  system 


m 

x  =  f{x)  +  Y  min  Hi.  u  =  (itijj  •  •  •  Um)  g  u 
2=1 

is  such  that  f  is  a  weakly  positively  Poisson  stable  vector  field,  and  accessibility  PARC 
is  satisfied,  then  the  system  with  controls  constrained  by  Uj  G  {  M;.  M;} .  M;  >  0,  i  = 

1,  •  •  • ,  m  is  controllable. 

In  the  setting  of  Lie  Poisson  reduced  dynamics  we  can  make  the  following  observation. 

Theorem  3  Let  G  be  a  Lie  group  that  acts  on  itself  by  left  (right)  translations.  Let 
H  :  T*G  — >■  1R  be  a  left  (right)  invariant  Hamiltonian.  Then, 

(i)  If  G  is  a  compact  group,  the  coadjoint  orbits  of  Q*  =  T*G/G  are  bounded  and  the  Lie- 
Poisson  reduced  Hamiltonian  vector  field  X~  defined  by  X~(ji)  =  {p,H}_(+j  is  WPPS. 

(u)  If  G  is  a  noncompact  group  then  the  Lie-P oisson  reduced  Hamiltonian  vector  field 
X~  is  WPPS  if  there  exists  a  function  V  :  Q*  — >■  R  such  that  V(p)  is  bounded  below, 
V(fi)  — »  oo  as  fi  — >■  oo  and  V  =  0  along  trajectories  of  the  system. 

Here  H  =  H  | g„  is  the  restriction  of  H  to  the  quotient  manifold  Q*  =  T*G/G  and  {•,  •}-(+) 
is  the  induced  minus  (plus)  Lie-Poisson  bracket  on  the  quotient  manifold  Q*  =  T*G/G. 

Proof:  (i)  The  map  A  :  T*G  — >■  Q*_  is  a  Poisson  map,  and  the  Poisson  manifold  Q*_  is 
symplectically  foliated  by  co-adjoint  orbits,  i.e.  it  is  a  disjoint  union  of  symplectic  leaves 
that  are  just  the  co-adjoint  orbits.  Any  Hamiltonian  system  on  Q*_  leaves  invariant  the 
symplectic  leaves  and  hence  restricts  to  a  canonical  Hamiltonian  system  on  a  leaf.  To 
study  the  dynamics  of  a  particular  system  with  initial  condition  //( 0)  G  Q*_,  we  therefore 
restrict  attention  to  the  co-adjoint  orbit  through  ^(0).  By  hypothesis,  each  co-adjoint 
orbit  is  compact.  The  flow  starting  at  ^(0)  preserves  the  symplectic  volume  measure  on 
the  orbit.  Hence  by  the  Poincare  Recurrence  Theorem,  we  know  that  for  almost  every 
point  p  G  Q*_  and  any  neighborhood  Vp  of  p  there  exists  a  time  t  >  T  such  that  (p) 
returns  to  Vp  i.e.  X~  is  WPPS. 

(ii) Let  D  =  {//  |  V(fi)  <  E}  and  Orb(-)  denote  the  coadjoint  through  p(0)  in  Q*_.  Then 
the  integral  curve  of  X~  starting  at  ^(0)  lies  entirely  in  the  set  S  =  iAnOrb(-).  Since  S 
closed  and  bounded  in  Q*_,  it  is  compact  in  Orb(-),  and  hence  as  before  X~  is  WPPS.  □ 

In  many  situations  the  function  H^,  =  H  +  (f>(Ci)  where  H  is  the  reduced  Hamiltonian 
and  Ci  are  the  Casimirs  are  a  good  choice  for  V  ( • ) . 

Remark:  In  our  present  setting  of  Lie-Poisson  reduced  dynamics,  WPPS  conditions  in 
theorem  2  can  be  verified  whenever  the  hypotheses  of  theorem  3  hold.  Once  WWPS  of  the 
drift  vector  field  has  been  established  theorem  2  can  be  used  to  conclude  controllability. 
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4  Examples 


In  this  section  we  discuss  the  controllability  of  the  Lie-Poisson  reduced  dynamics  of 
hovercraft,  spacecraft  and  the  underwater  vehicles.  These  systems  satisfy  conditions  of 
theorem  3.  The  kinematics  and  dynamics  of  these  examples  can  also  be  found  in  [1,  25,  26] 
for  completeness  we  present  necessary  details  here. 


4.1  Hovercraft  -  Planar  Rigid  Body  with  a  Thruster 

In  this  section  we  discuss  the  dynamics  of  a  planar  rigid  body  with  a  thruster.  The 
configuration  of  the  system  is  shown  in  Figure  1.  Let  { e\ .  .  ej  }  be  an  inertial  frame  of 

reference  fixed  at  O  and  {ej,  e|,  e| }  be  a  body  frame  fixed  on  the  rigid  body  B  at  its  center 
of  mass.  Since  the  rigid  body  is  restricted  to  move  in  the  e’jel]  plane,  a  typical  material 
point  qb  in  the  rigid  body  is  then  represented  in  the  inertial  frame  as  qr  =  Rqb  +  r  where 
R  is  an  element  of  SO( 2),  the  special  orthogonal  group  of  2  x  2  matrices  and  r  =  {x,y) 
is  a  vector  from  O  to  the  center  of  mass  of  B.  Hence  at  any  instant,  the  configuration 
X (t)  of  B  can  be  uniquely  identified  by  the  pair  (R,r)  or  equivalently  as  an  element  of 
SE( 2),  the  Special  Euclidean  group  of  2  x  2  matrices.  Recall 

SE( 2)  =  {(q  ’'I  |fteSO(2),re]R2} 

Let  us  assume  that  the  thruster  is  mounted  at  the  point  C  defined  by  the  vector  clb  in 
body  coordinates  and  dr  in  the  inertial  frame  of  reference.  The  thrusters  exert  a  force 
fr  in  inertial  coordinates  such  that  the  line  of  action  of  the  force  passes  through  C  and 
makes  an  angle  0  with  the  vector  clb.  We  now  derive  the  equations  of  motion  of  a  rigid 
body  subject  to  a  force  fr  along  a  specified  line  of  action. 


4.1.1  Symmetry  and  Reduction 

We  assume  for  now  that  the  rigid  body  (which  will  later  be  approximated  to  a  hovercraft) 
has  sufficient  lift  and  can  glide  on  the  surface  with  no  friction.  The  Lagrangian  L  : 
TSE( 2)  — >■  1R  for  this  case  is  simply  the  kinetic  energy,  i.e. 

L(i?,r,B,r)  =  i/n2  +  ^||r||2  (9) 

where  m  is  the  total  mass,  I  is  the  moment  of  inertia  of  B  in  the  body  frame,  is  the 
scalar  body  angular  velocity  about  the  center  of  mass.  The  corresponding  Hamiltonian 
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Figure  1:  Planar  rigid  body  with  thruster 


is  given  by 


H 


iPir 

2  m 


(10) 


where  II  =  JO  is  the  body  angular  momentum  and  p  =  mr  is  the  spatial  linear  momen¬ 
tum. 

Collecting  together  the  Newton-Euler  balance  laws  one  can  write  the  dynamics  in  spatial 
(inertial)  variables  (R,  r,  ir,p)  as 


i? 

=  M) 

(ll-a) 

r 

=  p/m 

(11-b) 

7T 

II 

a. 

X 

(11-c) 

P 

=  r 

(11-d) 

where 


We  observe  that  the  lifted  action  of  G  on  TSE( 2)  defined  by 

<f>9*  :  TSE{ 2)  — >■  TSE(2) 

( R.  r,  R.  f)  e- )►  (RR,  Rr  +  f,  RR,  Rr ) 
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leaves  the  Lagrangian  (9)  invariant.  Hence  the  Hamiltonian  is  also  G'- invariant.  We  can 
now  induce  a  Hamiltonian  on  the  quotient  space,  T* SE(2) / SE(2),  and  express  the  dy¬ 
namics  in  terms  of  the  appropriate  reduced  variables.  The  quotient  space  T* SE(2) / SE(2) 
is  isomorphic  to  se( 2)*,  the  dual  of  the  Lie  algebra  of  SE( 2)  and  the  reduced  variables 
are: 


n,  the  body  angular  momentum, 

P  =  RTp ,  the  converted  linear  momentum. 


The  reduced  Hamiltonian  H  is  given  by 


H=2ln  + 


2  m 


Choosing 


AG  = 


A,  = 


0 

0 

1  ^ 

\ 

( 

0 

0 

°\ 

0 

0 

0 

to 

II 

0 

0 

1  and 

0 

0 

0. 

/ 

V 

0 

0 

oy 

(12) 


as  a  basis  for  .se( 2)  we  have  the  commutation  relations:  [AG,  AG]  =  0,  [AG,  AG]  =  —AG  and 
[AG,  AG]  =  AG-  The  Lie-Poisson  bracket  of  two  differentiable  functions  G.  H  on  se(2)*  is 
then  given  by 

{G,P}_(p,)  =  X7GtA(/i)X7H  (13) 


where  //  =  (Pi,P2,n)  G  .se(2)*  and 


A 


0  0  P2 

0  0  -Pi 

-Pi  Pi  0 


The  reduced  equations  take  the  form 

i1,  —  {  /  /  ,  P  }  —  T  fext 


(14) 


where  fext  is  the  external  force  projected  appropriately.  In  the  present  setting  fext  = 
(|P|  cosG,  |P|  sinG,  |P||d|  sinG)T,  where  F  =  RT fr  and  |  •  |  denotes  the  norm  (see  [25]  for 
more  details). 
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Depending  on  the  control  authority  we  distinguish  two  versions  of  the  problem. 

Case  1:  The  Jet-Puck  Problem:  Here  we  assume  that  the  line  of  action  of  the  force 
is  fixed  (i.e.  (p  is  fixed)  but  its  direction  can  be  reversed.  Written  as  a  control  system, 
equations  (13)  take  the  form 

Pi  =  P2Ii/I  +  au 
P2  =  —  PiTl/I  +  j3u 
II  =  dJiu 

where  a  =  cos  (p,  (3  =  sin  <p  and  u  G  [1,-1]. 

Case  2:  The  Hovercraft  Problem:  Here  we  assume  that  we  now  have  control  over 
both  the  magnitude  of  the  thrust  and  <p.  The  equations  now  take  the  form 


Pi  =  P2n/I  +  Ui  cos(w2) 

(16-a) 

P2  =  -PiU/I+  ui  sin(w2) 

(16-b) 

n  =  dui  sin(w2) 

(16-c) 

where  11 1  G  [  1, 1]  and  u2  G  \pmin->  4*max\ 

(15-a) 

(15-b) 

(15-c) 


4.1.2  Controllability  of  the  Reduced  Dynamics 

Proposition  1:  The  jet-puck  dynamics  defined  by  (9)  are  controllable. 
Proof:  We  first  show  that  LARC  is  satisfied.  To  show  that 


dim(span£{yiS})(p)  =  3,  Vp  €  se(2)* 
where  /  =  (P2U/I,  —  P\Ii/I,  0)T  and  g  =  (o,  /3,  cl/3)T,  observe  that 


det {g,  [[/, g],g],  [[/, g],  [[/, g],g]\) 
a  2 7  /'  —2j j/32a 

=  det  I  *  d 


(3  —2  j(3a 

L  d/3  0 

o2) 


0 


=  -4^(;d2 

I  (d/?)4 


=  — 4- 


/ 3 


(since  ft2  +  / 3 2  =  1) 


Hence  dim(span£{j!g|)(p)  =  3  Vp  G  se(2)*  as  long  as  / 3  =  sin  ip  ^  0,  i.e.  as  long  as  the 
line  of  action  of  F  does  not  pass  through  the  center  of  mass. 
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Figure  2:  Energy  surface  and  coadjoint  orbits  in  se(2)* 


We  observe  that  the  reduced  Hamiltonian 


H 


2 


2  m 


(17) 


is  bounded  below,  radially  unbounded  and  is  such  that  H  =  0.  Hence  it  follows  from 
theorem  3  that  /  is  WPPS  and  hence  from  theorem  2  we  conclude  that  the  jet-puck 
dynamics  are  controllable.  □ 

In  fact  we  one  observes  that  every  orbit  of  /  is  periodic  and  hence  trivially  Poisson  stable. 
Remark:  Observe  that  the  coadjoint  orbits  in  se(2)*  are  cylinders 

{(Pi,  P2,  n)  6  5R3  I  P\  +  P22  =  constant  ^  0}. 

_  p2  p2  yj2 

The  surfaces  defined  by  D  =  {(Pi,P2,n)  |  2m  +  2m  +  =  const}  are  ellipsoids.  From 

theorem  3  the  integral  curves  of  the  the  vector  field  are  restricted  to  the 

set  S  =  PnOrb(-),  which  in  this  case  is  simply  S 1  (see  fig  (2)). 

Proposition  2:  The  hovercraft  dynamics  defined  by  (10)  are  controllable. 

Proof:  In  (10)  setting  u2  =  A:,  where  A:  is  some  constant  not  equal  to  zero,  the  equa¬ 
tions  reduce  to  those  of  the  jet-puck  and  hence  from  Proposition  1  the  dynamics  are 
controllable.  □ 
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4.1.3  Cotangent  Space  Controllability 


We  now  investigate  the  controllability  in  the  cotangent  space.  We  show  that  the  complete 
(unreduced)  system  is  locally  strongly  accessible.  The  complete  set  of  equations  for  the 
jet-puck  take  the  form 


X 

=  pi /m 

(18-a) 

y 

=  p2/m 

(18-b) 

6 

=  n// 

(18-c) 

Pi 

=  (cos  (0  +  (f)))u 

(18-d) 

P2 

=  (sin(0  +  <p)u 

(18-e) 

ri 

=  (d  sin  cj))u 

(18- f) 

Proposition  3  :  For  the  system  defined  by  (17)  the  cotangent  space  T*SE( 2)  is  locally 
strongly  accessible. 

Proof:  To  show  that 

dim(span£{j!9})(p)  =6,  Mp  G  T*SE( 2) 
where  /  =  (pi/m,p2/m,  II//,  0,  0,  0)T, 

and  g  =  (0,  0,  0,  cos (9  +  </>),  sin(#  +  </>),  cl  sin  </>)T,  observe  that 

=  16/v  (cos2 (6  +  </>)  +  sin2(6*  +  </>))2 
=  16  K 


where 


6  =  [f,g],  6  =  [[f,g],g],  6  =  [[[f,g],f],g], 
Ci  =  [[[[f,g],f],g],  [/w]], 

&  =  [[[[/,  0], /],$],  [[/w]w]]  and 

W  =  (i)2(I)7(dsin^)8 


Hence  again  if  sin^  ^  0,  dim(span£{ji9})(p)  =  6,  \/p  G  T*SE( 2).  Also  [/,  JY]  G 

span((7,  ^1,  £2,  £3,  £4,  £5),  V-Y  G  {g,  6, 6, 64, 65}-  Hence  the  complete  system  is  locally 

strongly  accessible.  □ 
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4.1.4  Small-Time  local  controllability 


Recall  that  a  control  system  (8)  is  said  to  be  small  time  locally  controllable  (STLC)  from 
x0  £  M  if  it  is  locally  accessible  from  x0  and  x0  is  in  the  interior  of  Rl  (x0,  <  T)  1  for  all 
T  >  0  and  each  neighborhood  V  of  x0.  If  this  holds  for  any  x0  £  M  then  the  system  is 
called  STLC. 

In  [27]  Sussmann  proves  the  following  sufficient  condition  for  scalar  input  systems. 
Theorem  4  Consider  an  analytic  system 

x  =  fo{x)  +  fi{x)u  [u{t)\  <  A  (19) 

and  a  point  x0  such  that 


where  S1(Xi,X2)  is  the  linear  span  of  X1,X2,  and  the  brackets  (adXi)C X2  for  j  >  1  and 
{u  |  fo{xo )  +  ufi(xo)  =  0}.  Then  (19)  is  not  STLC  from  xq- 

We  use  the  above  result  to  show  that  the  jet-puck  dynamics  are  not  STLC 

Proposition  4  The  jet  puck  dynamics  defined  (17)  are  not  STLC  from  the  origin. 

Proof:  It  is  sufficient  to  consider  STLC  of  the  reduced  dynamics.  With  u  =  0  ob¬ 
serve  that  S1(f,g)( 0)  is  a  one-dimensional  space  spanned  by  while 

\9,  L/w]](0)  =  -2jf32^  +  Hence  [g,  [/,  #]](0)  ^  S1{f,  ^)(0)  □ 

4.2  Attitude  control  of  spacecraft 

We  now  discuss  the  dynamics  describing  spacecraft  attitude  control  with  gas  jet  actuators. 
Let  be  a  body  frame  fixed  on  the  rigid  body  (spacecraft)  B  at  its  center  of 

mass  and  let  { cj .  ef .  of  }  be  an  inertial  frame  of  reference  with  origin  coincident  with  the 
origin  of  the  body  fixed  frame  (see  Fig  3).  A  typical  material  point  qb  in  the  rigid  body 
is  then  represented  in  the  inertial  frame  as  qr  =  Rqb  where  R  is  an  element  of  SO(3 ),  the 
special  orthogonal  group  of  3  x  3  matrices.  Hence  the  configuration  space  of  the  rigid 
body  may  be  identified  with  an  element  of  SO(3 ),  the  velocity  space  is  with  the  tangent 
bundle  TSO( 3)  and  the  momentum  phase  space  with  the  cotangent  bundle  T*SO( 3). 
Let  b\,---bm  be  the  axis  about  which  the  corresponding  control  torque  of  magnitude 

1Rl  (xo,  <  T)  =  |J  Rl  (xo ,t)  where  Rl  ( xq,T )  is  the  reachable  set  from  .To  at  time  T  >  0,  following 

0  <t<T 

trajectories  which  remain  for  t  <  T  in  the  neighborhood  V  of  tq 
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Figure  3:  Spacecraft  with  gas  jets 

|  b;  j  //,  is  applied  by  means  of  opposing  pairs  of  gas  jets.  The  dynamic  equations  for  the 
controlled  space  craft  are  then  given  by 

R  =  RQ  (20-a) 

m 

IQ  =  IOxO  +  ^Mi  (20-b) 

2  =  1 

where  Q  =  (fix,  Q2.  fi3)  is  the  body  angular  velocity,  Q  is  a  3  x  3  skew  symmetric  matrices 
given  by 

o  -n3  n2 

n  =  o3  o  -fii 

—n2  n  i  o 

and  /  =  diag(/| .  J2,  /3)  is  the  inertia  matrix.  In  the  rest  of  the  discussion  '  defines  a  map 
1R3  — >■  so( 3),  such  that  S/3  =  a  x  /3,  a ,  0  G  1R3.  Thus 

0  — ct3  ol2 

ex  =  a3  0  — tti 

—  tti  0 

4.2.1  Symmetry  and  Reduction 

The  Lagrangian  L  :  TSO( 3)  — >•  H»  is  again  simply  the  kinetic  energy  and  is  given  by 

l{r,r)  =  -  <  njn  > 

2 
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and  the  corresponding  Hamiltonian  H  :  T*SO( 3)  — >■  1R  is  given  by 

\  <  n,/_1n  > 

where  n  =  /fi  is  the  body  angular  momentum.  Observe  that  the  lifted  action  of  g  = 
R  G  50(3)  on  T50( 3)  defined  as 

:  TSO{3)  ->•  T50(3) 

(/?.  i?)  I—?-  (iLR,  RRQ) 

the  lifted  action  of  g  =  .R  G  50(3)  on  TSO(3)  defined  as  leaves  the  Lagrangian  (and 
hence  also  the  Hamiltonian)  invariant.  Hence  one  can  induce  a  Hamiltonian  on  the 
quotient  space,  T*  SO  (3)  /  SO  (3),  and  express  the  dynamics  in  terms  of  the  appropriate 
reduced  variables.  The  quotient  space  T*50(3)/50(3)  is  isomorphic  to  so( 3)*,  the  dual 
of  the  Lie  algebra  of  50(3)  and  the  reduced  variables  are  n  =  (Hi,n2,n3)  the  body 
angular  momentum.  Choosing 


1 

( 

0 

0 

0 

\ 

(° 

0 

1 

II 

0 

0 

-1 

to 

II 

0 

0 

0 

1 

K 

0 

1 

0 

) 

1 

l-l 

0 

0 

( 

0 

— ; 

1  0 

\ 

Ad  = 

1 

0 

0 

\ 

0 

0 

0 

) 

and 


as  a  basis  of  .so(3)  and  with  the  commutation  relations  [Ad,  A"2]  =  X3 ,  [A3,  Ad]  =  Ad  and 
[Ad,  Ad]  =  —Ad,  the  Lie-Poisson  bracket  of  two  differential  functions  O.  //  on  so(3)*  is 
given  by 


(G,H}_(/u)  =  VGtA(/u)VR 


where  //  =  (Hi,n2,n3)  G  so(3)*  and 


A 


o  —  n3  n2 
n3  o  —Hi 
— n2  n3  o 


The  reduced  equations  take  the  from 


.  hh  J  m 

ri2  =  l^u2u3  +  h'u' 
ri3  =  hr^nlu2  -1 


(21) 


(22) 
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where  6* 

=  RTb;r  Assuming  that  we  have  only  one  control  and  I\ 

=  J2  (21)  can  be  written 

as 

rii 

=  (/l~/3)n2n  3  + cm 

iiH 

(23-a) 

ri2 

=  -{h-h)ThTh  +  f3u 

4 1^3 

(23-b) 

ri3 

=  7  u 

(23-c) 

4.2.2  Controllability  of  reduced  dynamics 

Observe  the  similar  structure  of  base  space  equations  for  the  jet-puck  and  those  of  the 
controlled  Euler  equations.  Hence  we  can  make  similar  claims  regarding  controllability 
and  STLC.  Proofs  are  omitted  as  they  are  similar  to  those  of  the  jet-puck  dynamics. 

Proposition  5  The  spacecraft  dynamics  defined  by  (22)  are  controllable  if  a 2  +  f32  ^  0 
and  7  ^  0 

Proposition  6  The  spacecraft  dynamics  (22)  are  not  STLC  from  the  origin. 

Remark:  The  coadjoint  orbits  in  so(3)*  are  spheres  (see  fig(4)) 

{(nx, n2, n3)  e  r3  |  +  ng  =  const.}. 

In  this  case  since  the  coadjoint  orbits  are  compact  manifolds  one  can  conclude  from 
theorem  2  that  the  drift  vector  field  is  WPPS.  Fig  (4)  shows  the  intersection  of  the 
coadjoint  orbits  and  the  energy  surface. 


4.3  Autonomous  Underwater  Vehicle 

In  this  section  we  discuss  the  reduced  space  controllability  for  a  neutrally  buoyant  un¬ 
derwater  vehicle  (UV).  We  distinguish  between  the  cases  of  coincident  and  noncoincident 
centers  of  buoyancy  and  gravity.  The  Lie-Poisson  dynamics  for  these  cases  have  been 
derived  in  [26].  We  only  present  a  brief  overview  of  the  Lie-Poisson  dynamics  and  then 
use  the  results  of  theorem  2  and  theorem  3  to  make  necessary  conclusions  about  control¬ 
lability. 

Let  {cj'.cT.  c3}  be  an  inertial  frame  of  reference  (see  fig  5)  fixed  at  O  and  {e5,e2,e|}  be 
a  body  frame  fixed  on  the  vehicle  at  its  center  of  buoyancy  (CB).  A  material  point  qb  in 
the  UV  is  then  represented  in  the  inertial  frame  as  qr  =  Rqb  +  r  where  R  is  an  element  of 
SO( 3),  the  special  orthogonal  group  of  3  x  3  matrices  and  r  =  (x,y,z)  is  a  vector  from 
O  to  the  center  of  buoyancy  (CB).  Hence  at  any  instant,  the  configuration  X(t)  of  the 
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b 


Figure  4:  Energy  surface  and  coadjoint  orbits  in  se( 2)* 


UV  can  be  uniquely  identified  by  the  pair  ( R,r )  or  equivalently  as  an  element  of  SE( 3), 
the  Special  Euclidean  group  of  3  x  3  matrices.  Recall 

SE{  3)  =  {^  l^J  I  R  e  SO(3),r  e  IR3} 

While  deriving  the  dynamics  we  assume  that  the  UV  is  submerged  in  an  infinitely  large 
mass  of  incompressible,  inviscid  fluid.  Further,  we  assume  that  the  flow  is  irrotational 
(the  motion  of  the  fluid  is  entirely  due  to  that  of  the  UV).  Under  these  assumptions  the 
motion  of  the  fluid  can  be  characterized  by  the  existence  of  a  single  valued  potential  <f> 
which  satisfies 

V2</>  =  0 
VR  =  0  at  infinity 

dcj) 

—  — —  =  n.(v  +  0  x  rb)  at  body  surface, 


where  rb  is  a  vector  from  the  CB  to  the  vehicle’s  surface,  n  is  the  unit  outward  normal  vec¬ 
tor  of  the  vehicle,  =  (IR,  U2,  f23)T  are  the  body  angular  velocities,  and  v  =  (rq,  vz)T 
are  the  linear  velocity  components  along  the  body-fixed  frame.  Under  these  assumptions 
Kirchhoff  showed  that 


<j)  =  Vl(j>i  +  v2h  +  v3(j)3  +  Q.1X1  +  ^2X2  +  ^3X3  (24) 

where  Ri,  </>2,  <j> 3,  Xi)  X2,  X3  are  functions  of  x,  y,  z  determined  by  the  configuration  of  the 
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Figure  5:  Autonomous  Underwater  Vehicle 


surface  of  the  solid.  Using  the  form  of  d>  as  expressed  in  (23)  the  kinetic  energy  of  the 
fluid 


where  p0  is  the  fluid  density,  can  be  expressed  as  a  quadratic  from  as 

Tf  =  5H’T0H'  0  =  (  Z 


©12  \ 
©22  J 


0ii  is  referred  to  as  the  added  mass  matrix  and  ©22  as  the  added  inertia  matrix 

Assume  that  the  center  of  gravity  (CG)  does  not  coincide  with  the  center  of  buoyancy 
(CB)  and  lies  on  the  e3  axis  at  a  distance  l  >  0  (bottom  heavy)  from  the  CB  i.e  rg  =  U3 
where  i 3  denotes  a  unit  vector  (in  body  coordinates)  along  the  e3  axis.  Also  let  ig  denote 
a  unit  vector  (in  inertial  coordinates)  in  the  direction  of  gravity,  i.e.  along  the  e3  axis. 
Let  m  be  the  mass  of  the  vehicle,  Jb  the  inertia  matrix  for  the  vehicle.  The  Lagrangian 
L  :  TSE(3)  — >■  R  is  then  given  by 

L(R ,  r,  R ,  r)  =  -(0T  JQ  +  2 QT Dv  +  vT Mv  +  2mgl{ig.Ri^)) 

Z, 

where  J  =  Jb  +  0n,  D  =  mlfg  +  0i2  and  M  =  ml  +  022  (I  is  the  3  x  3  identity 
matrix).  In  the  rest  of  the  discussion  the  UV  is  approximated  as  an  ellipsoid  and  hence 
©12  =  ©21  =  0. 
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Observe  that  the  Lagrangian  is  invariant  under  the  action  of  the  group 

G  =  {(if,  r)  e  SE{ 3)  |  RTig  =  ig}  =  SE{ 2)  x  R. 

and  hence  the  Hamiltonian  system  on  T*SE( 3)  (which  is  also  left  invariant  under  the 
action  of  SE( 2)  x  R)  can  be  reduced  to  a  Hamiltonian  system  on  <S*,  the  dual  of  the  Lie 
algebra  of  the  semi-direct  product  S  =  SE( 3)  xp  R3  (see  [26]  for  details).  The  reduced 
Hamiltonian  on  S*  is 

H{ n,  p,  , )  =  i(nTon  +  2nTpTp  +  ptcp  -  2 mgi{,  ,?;3)), 

z 


where 


A  =  (J  —  DM~1Dt)-\  B  =  -CDTJ~\  C  =  (. M  -  DtJ~1D)-\ 
n  =  .70  +  Dr.  P  =  Mv  +  DT O,  and  ,  =  RTig. 


Choosing 


Bi 


A;  0\ 

0  0  )  • 


1,  -  -  - ,  6,  B., 


0  &i—6 

0  0 


?'  =  7,  8,  9. 


where 


A  = 


0 

0 


i  =  1,2,3  A 


0  e,; 

0  0 


/  =  4,  5,  6 


as  a  basis  for  <S  the  Lie  algebra  of  S  the  Lie-Poisson  bracket  of  two  differentiable  functions 
G,  77  on  S*  is  given  by 

{<A77}_(A  =  VGtA(AV77 


where  /.i  =  (n,  P, ,  )  and 


A 


n  p 
P  0 


0 

0 


The  Lie-Poisson  reduced  equations  (see  [26]  for  a  complete  description  of  reduction  pro¬ 
cedure)  are  then  given  by 


Ik  =  {Ah,P}-(fO 
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or  explicitly  as 


n  x  (An  +  BtP)  +  P  x  (CP  +  pn)  —  mgl ,  x  i3 

(25-a) 

p  x  (An  +  btP) 

(25-b) 

,  x  (An  +  btp) 

(25-c) 

Proposition:  The  Lie-Poisson  reduced  Hamiltonian  vector  field  (given  by  (25))  defined 
on  5*  is  WPPS. 

Proof:  Choose  T  (n,  P, ,  )  =  H( n,  P, ,  )  +  ,  T,  .  Observing  that  V  is  radially  unbounded 
and  that  V  =  0  along  trajectories  of  (25)  the  result  follows  form  theorem  3.  □ 

Remark  :  In  the  case  of  coincident  center  of  gravity  and  center  of  buoyancy  (i.e.  1  =  0), 
the  Hamiltonian  system  on  T*SE( 3)  is  left  invariant  under  the  SE( 3)  action  of  rotations 
and  translations,  and  we  can  derive  a  set  of  reduced  Lie-Poisson  equations  on  se(3)*. 
Choosing 


A  = 


o 

0 


i  =  1,2,3  /  A 


0  e'i 

0  0 


i  =  4,  5,  6 


as  the  basis  for  .se(2)  the  structure  matrix  [A(/i)]?y  =  —  Y.l=icij^k  is  given  by 


a(aO  =  A(n,  P) 


n  P 
P  o 


where  n  =  JO  and  P  =  Mv.  The  Lie-Poisson  reduced  equations  are  given  by 

(26-a) 
(26-b) 


ri  =  n  x  (An)  +  P  x  CP 

P  =  Px  An 


Choosing  V(n,P)  =  H  =  |(nTAn  +  PTCP),  the  reduced  Hamiltonian,  we  observe  that 
the  Lie-Poisson  reduced  Hamiltonian  vector  field  defined  on  se(3)*  is  WPPS. 


4.3.1  Controllability  of  Reduced  Dynamics 


Assuming  that  we  have  three  controls  U\.  u2,  «3  such  that  U\  provides  a  pure  torque 
about  c, .  U2  provides  a  pure  torque  about  e|  and  u?t  provides  a  pure  translation  along 
the  e\  axis.  The  Lie-Poisson  reduced  dynamics  with  controls  for  the  underwater  vehicle 
with  coincident  CB  and  CG  are 


A 


h-h 

hh 


n2n3  + 


m2  —  m  3 
m2m3 


P2P3  +  U  | 


(27-a) 
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ri2 

ri3 

A 

P2 


h~h 
hh 
h  ~  1 2 

hh 

p2n3 


narix  + 

rprp  + 

-P3  n2 


h 

1 

1 - 1 

to 

P3  ip 

Pi  n3 

h 

h 

P1P2 

P2n, 

1 — >H 

to 

h 

777-  3  —  777-1 
777-3  777-1 
777-1  —  ?n2 
777-1 777-2 

I-  «3 


P3P1  +  77-2 

PlP2 


(27-b) 

(27-c) 

(27-d) 

(27-e) 

(27-f) 

(27-g) 


Proposition:  The  Lie-Poisson  reduced  dynamics,  defined  by  (27),  of  the  underwater 
vehicle  with  coincident  center  of  buoyancy  and  center  of  gravity  are  controllable  if  I\  ^  h- 

Proof:  Observing  that 


deb{gugi,gz,  [[/, £1], g2\,  [[/,  #2],  £3],  [[[/,02],  [Jws]],  £1])  =  ^2^  P  0 


if  I\  7^  P  where 


/  = 


P|P2 


/  Ifc^iprLs  +  mp  p 

/  /2i3  z  0  m2m3  z  0 

iaz_hn3ni  +  p 

I3I1  °  x  m3mi  0  x 

+  2*1=2* 

J1J2  toito-2 

P2  III  _  P3II2 

7s  h 

P31I1  _  Pin3 

71  73 
Pi  n2  P2  ft  1 

72  7i 


( 1 ) 

(  0  ^ 

(  0  \ 

0 

1 

0 

0 

0 

0 

,  Ui  = 

0 

,  u2  = 

0 

,  U\  = 

1 

0 

0 

0 

V  0  j 

V  0  ) 

1 0  / 

i.e.  dim(span£{yiflli92)g3})(p)  =  6,  \/p  €  se(3)*,  and  that  /  is  WPPS  the  result  follows 
from  theorem  2.  □ 


5  Conclusions  and  Future  Work 


We  have  discussed  the  controllability  of  the  Lie-Poisson  reduced  dynamics  of  a  class  of 
mechanical  systems  and  present  an  approach  to  conclude  controllability.  The  Lie-Poisson 
reduced  dynamics  and  controllability  of  hovercraft,  spacecraft  and  underwater  vehicle  are 
discussed. 

Current  research  includes  the  design  of  feedback  laws  to  stabilize  the  origin  of  the  reduced 
system.  We  are  also  studying  the  use  of  periodic  controls  to  generate  loops  in  the  base 
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space,  and  thereby  steer  in  the  fiber.  Some  future  research  includes  control  laws  to  steer 
hovercraft  and  hybrid  architectures  for  generation  of  “scripts”  for  obstacle  avoidance  and 
navigation. 
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Appendix 

A  Accessibility  and  Controllability 

We  briefly  review  definitions  and  related  theorems  on  accessibility  and  controllability  of 
smooth  affine  nonlinear  systems, 


x  =  f{x)  +  Y,g^x)ui,  (28) 

2  =  1 

where  x  =  (.To,  •  •  • ,  x„  )  are  local  coordinates  for  a  smooth  manifold  M  and  u  =  (u\,  ■  ■  ■  um)  G 
U  C  !Rm.  We  follow  the  development  in  [19]  (see  historical  references  therein).  It  is  as¬ 
sumed  that  - 

(i)  The  input  space  U  is  such  that  the  set  of  associated  vector  fields  of  (28)  T  =  {/( x)  + 
Y.T=i  //,(./')//;  ( // 1 .  •  •  •  Um)  G  U }  contains  the  vector  fields  /,  gi,  ■  ■  ■ ,  gm. 

(ii)  The  set  of  admissible  controls  consists  of  piecewise  constant  functions  which  are 
piecewise  continuous  from  the  right. 

Let  Rl  (t0,T)  be  the  reachable  set  from  x0  at  time  T  >  0,  following  trajectories  which 
remain  for  t  <  T  in  the  neighborhood  V  of  x0.  Let  Rl  (t0,  <  T)  =  |J  Rl  (x 0,t) 

0  <t<T 

Definition:  The  system  (28)  is  locally  accessible  from  t0  if  Rv  (t0,  <  T )  contains  a 
non-empty  open  set  of  M  for  all  neighborhoods  V  of  x0  and  all  T  >  0.  If  this  holds  for 
any  To  G  M  the  the  system  is  called  locally  accessible. 

Let  L(x)  =  {spamY(.T)  |  X  vector  field  in  £,  x  G  M)  where  £,  the  accessibility  Lie 
algebra,  is  the  smallest  subalgebra  of  the  Lie  algebra  of  vector  fields  on  M  that  contains 
/,  <7i,  •  gm-  The  system  is  said  to  satisfy  the  accessibility  Lie  algebra  rank  condition 
(LARC)  if 

L(x)  =  TXM ,  Vt  G  M  (29) 

Theorem  5  If  dimL(x)  =  n  Vt  G  M  then  the  system  system  (28)  is  locally  accessible. 

Definition:  The  system  (28)  is  said  to  be  locally  strongly  accessible  from  t0  if  for  any 
neighborhood  V  of  To  the  set  Rl  (xq,T)  contains  a  non-empty  open  set  for  any  T  >  0 
sufficiently  small. 

Let  C0  be  the  smallest  Lie  subalgebra  which  contains  g\ ,  •  •  • ,  gm.  and  satisfies  [f,X]  G 
Co.  VA  G  £o  and  L0(x)  =  {span A (t)  |  X  vector  field  in  Cqs  x  G  M}.  The  system  is  said 
to  satisfy  the  strong  accessibility  Lie  algebra  rank  condition  if 

L0(x)  =  TXM ,  Vt  G  M  (30) 
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Theorem  6  If  dimL0(x0 )  =  n,  then  the  system  (28)  is  locally  strongly  accessible  from 
x0. 


Let  x(t,0,  x0,u)  denote  the  solution  of  (28)  at  time  t  >  0  for  a  particular  input  function 
u ( • )  and  initial  condition  rc(0)  =  x0- 

Definition:  The  system  (28)  is  called  controllable  if  for  any  two  points  xi,x2  in  M  there 
exists  a  finite  time  T  and  an  admissible  function  u  :  [0,T]  — >■  U  such  that  x(t,0,xi,u)  = 
x2. 


Theorem  7  (Chow) [28]  The  nonlinear  system 


m 

x  =  9i(x)Ui,  U  =  {uU  ■  ■  ■  ,  Um)  C  U 
i— 1 


is  co7itrollable  if  the  accessibility  LARC  is  satisfied. 


(31) 
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